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ABSTRACT 

We investigate the evaporation history of known transiting exoplanets in order to 
consider the origin of observed correlations between mass, surface gravity and orbital 
period. We show that the survival of the known planets at their current separations 
is consistent with a simple model of evaporation, but that many of the same planets 
would not have survived closer to their host stars. These putative closer-in systems 
represent a lost population that could account for the observed correlations. We con- 



clude that the re lation underlying the correlations noted by Mazeh et al. (20051 and 



Southworth et al.| ( |2007[ ) is most likely a linear cut-off in the M'^/R^ vs plane, 
and we show that the distribution of exoplanets in this plane is in close agreement 
with the evaporation model. 
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1 INTRODUCTION 



When Mayor & Queloz ( 1995 I discovered the first planet 



around a main sequence star, its small orbital separation 
came as a surprise to the astronomical community, who had 
expected to find planets similar to those in our solar sys- 
tem. This motivated work into how and where these 'hot 
Jupiters' form and how they might migrate to their current 
separations (see for example Lin et al.|1996 Boss|1996 |. 

Many close-in planets are now known and intriguingly 
some planet properties seem to be correlated. |Mazeh et al.| 
( 2005 1 suggest that planet mass is correlated with orbital 



period, while Southworth et al. ( 2007 1 find that surface grav- 



ity may be correlated with orbital period. These correlations 
have been re-plotted by various authors as more planets have 



been discovered (for example Hansen &: Barman||2007 Pol- 
lacco et al.|2008 l. Both distributions are reproduced in our 



Figure [Ij using the sample of transiting planets from this 
work (see Section 



2.11 



It has been suggested that these correlations might be a 
result of the migration mechanism (e.g. [Torres et al.| ( |2008[ ) 
and references therein) or that they are a result of evapora- 
tion of close-in planets, such as that reported to be ongoing 
in HD209458b by[VidaI-Madjar et aL| ( [2003| (see also[Ben^ 
Jaffel|2007[|Vidal-Madjar et al.|2008| 



Lecavelier Des Etangs ( 2007 1 developed a simple model 



of exoplanet evaporation driven by X-ray/EUV radiation 
and considered the current evaporation rates of known 
planets. He found that known planets will not be losing 



significant mass at today's irradiation rates. |Penz et al.| 
(2008| ) pointed out that X-ray/EUV irradiation will be much 
stronger around young stars and showed that, in principle, 
the mass distribution of close-in exoplanets could be signif- 
icantly modified by evaporation. 

In this paper we use a simple model of evaporation to 
investigate whether the trends in planet properties identified 
bylMazeh et al.l (120051) andlSouthworth et al.|(|2007l) can be 



accounted for by past evaporation. We find good agreement 
between the distribution of known planets and the cut-off 
predicted by the evaporation model, and we take this as 
evidence that a population of close-in exoplanets has been 
lost to evaporation. 



2 METHOD 

2.1 Selection of exoplanet sample 

We select only transiting planets for our study because we 
require the planetary radius in order to consider surface 
gravity and evaporation rates. We selected planets from The 
Exoplanets Encyclopedicj^on 15 June 2008. Key parameters 
were selected from the encyclopedia and, where necessary, 
original discovery papers. These parameters allow the calcu- 
lation of the energy incident upon the planet per unit time, 
and the total potential energy of the planet. Stellar quanti- 
ties are shown in Table [l] and planetary properties in Table 
[2] for the sample of planets used in this work. 

We used as many of the known transiting systems as 
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Figure 1. Plots of the property correlations of |Ma2eh et al.| | [2005[ l 
and [Southworth et aT] | |2007[ l , updated to include the entire sam- 
ple of planets as listed in Table |2] Planets exist beyond the area 
of the plot, specifically HAT-P-2b, but always towards the top 
right. 



possible, although some were removed due to the absence 
of data required in our analysis. The CoRoT systems and 
Lupus- TR-3b were removed due to lack of distance infor- 
mation available at the time of our analysis. OGLE-TR- 
10b, 111b, 182b, 211b, and SWEEPS-04b and lib were re- 
moved due to the host stars being of unknown spectral type. 
GJ 436b was not included because it has the only M type 
host, and will have different high energy emisison character- 
istics to the FGK hosts prevalent in the rest of the sample. 
HD149026b was not included because its host star is not on 
the main sequence. 

In Figu re [l] we plot m ass against orbital period for our 
sample af ter [Mazeh et ah] ([2 005), and surface gravity vs pe- 
riod after Southworth et aL| ( |2007| ). Both plots show greater 
scatter than the original suggested correlations, especially 
to the upper right, but the mass-period distribution still 
shows a sharp cut-off to the lower left. The lower left of the 
surface gravity-period distribution does not show a sharp 
cutoff, but this area remains relatively sparsely populated, 
despite strong selection effects in transit surveys for plan- 
ets in this region. The upper right of both of these diagrams 
also remain sparsely populated. This is unlikely to be related 
to evaporation and instead may be a feature of planetary 
migration. Selection effects in transit surveys may also be 
a factor, since they select against long orbital periods and 
high surface gravities. 



It has been shown by Lammer et al. (20031 that Jeans es- 



cape of material from hot Jupiters cannot drive significant 
mass loss. Instead they suggest that energy-limited escape 
of particles due to absorption of EUV/X-ray photons is the 
most-likely driver of mass loss in these systems. 

In order to estimate the evaporation rates of exoplanets 



we adopt the model of Lecavelier Des Etangs (20071, who 



used the ratio of planetary potential energy and the time- 
averaged power incident on the planet in order to predict 
mass loss rates, given by 



3GMp{ay 



(1) 



where rh is mass loss rate, is the star's X-ray/EUV lumi- 
nosity, Rp is the planetary radius, Alp is the planetary mass, 
G is the gravitational constant, and (a) is the time-averaged 
orbital distance. 

This simple model assumes energy-limited evaporation, 
which may overestimate real mass loss rates. It also assumes 
that a test particle must be moved to infinity to escape the 



gravitational potential of a planet (Sect. 4.3 of Lecavelier 
[Pes Etangs, 2007p , which Erkaev et al. (20071 suggest may 
underestimate mass loss by up to 40% when a planet is close 
to filling its Roche lobe. Further, we neglect any evolution 
of the planet radius, assuming the current value for each 
planet. Since unperturbed gaseous planets contract slowly 
with time, this will tend to underestimate mass loss rates. 



2.3 Estimation of X-ray/EUV irradiation 

X-ray/EUV emission from stars is believed to be driven 
by a dynamo action that decreases with time as the star 
spins down due to magnetic braking ( Skumanich||1972 \ . X- 
ray/EUV emission is seen to increase with stellar rotation 
rate, up to a maximum of L^/Lboi — 10~^, where it satu- 
rates ( |Vilhu||1984[ [Vilhu fc Walter||1987[ ). We expect there- 
fore all known exoplanets to have suffered a period of high 
and approximately constant X-ray/EUV irradiation in the 
past, when their host stars were young and rapidly rotating. 
This saturated phase is likely to dominate their irradiation 
history. 

The duration of the saturated phase of X-ray emission 
depends on spectral type, but is not well known. For F type 
stars a saturation time of 0.1 Gyr is used, taken from |Flac^ 
comio et al. ( 2003 1 who used optical and X-ray telescope 



data of five star-forming regions of varying age (p Ophiuchi, 
the Orion Nebula Cluster, NGC 2264, Chamaeleon I, and 
rj Chamaeleontis) and two young clusters (the Pleiades and 
NGC 2516). This is their upper limit for the saturation pe- 
riod of stars with the mass of F stars. 

For G type stars we estimate a saturation time of 0.2 
Gyr using the results of Giidel et al. (19971, who used the 



ASCA and ROSAT X-ray satellites to probe the coronae of 
a sample of nine solar-like G stars, with ages between 70 
Myr and 9 Gyr. 

The saturation timescale for K type stars is estimated 
from the observation (e.g., Wheatley|1998 1 that K stars are 
saturated in the Pleiades (age of 0.1 Gyr, Basri et al.|1996 l 
but not in the Hyades (age of ~0.6 Gyr, [Lastennet et al.| 
20001. We choose a value of 0.35 Gyr, but adjusting this 
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Table 1. Stellar parameters used in this work. Stellar mass, radius and spectral type are taken from 
The Exoplanot Encyclopedia" unless otherwise noted, and bolometric magnitudes from [AUcn (20001. 
p, is calculated using the standard formula from the data listed. -LJ"* is calculated by converting the 
bolometric magnitude into a bolometric luminosity, and using the relations from |Pizzolato et al.| l [2003[ l 
as described in Section |2.3| Tsat is estimated as described in Section |2.3| 



Star Mass Radius p* Spec. Type Bol.Mag. logLJ"' Tsat ref 

(Mq) (Hq) gcm-3 (mag) (ergs s"!) (Gyr) 



HAT-P-1 


1.12 


1.15 


1.04 


GOV 


4.2 


29.9 


0.2 


a 


HAT-P-2 


1.298 


1.41 


0.65 


F8 


3.8 


30.06 


0.1 


a 


HAT-P-3 


0.936 


0.82 


2.36 


K 


5.6 


30.14 


0.35 


a 


HAT-P-4 


1.26 


1.59 


0.44 


F7V 


3.5 


30.18 


0.1 


a 


HAT-P-5 


1.16 


1.17 


1.03 


GIV 


4.5 


29.78 


0.1 


a 


HAT-P-6 


1.29 


1.46 


0.59 


FIV 


3.2 


30.3 


0.1 


a 


HD 17156 


1.2 


1.47 


0.53 


GO 


4.2 


29.9 


0.2 


a 


HD 189733 


0.8 


0.75 


2.64 


K1-K2 


5.6 


30.04 


0.35 


a 


HD 209458 


1.01 


1.12 


1.01 


GOV 


4.2 


30.6 


0.2 


a 


OGLE-TR-113 


0.78 


0.77 


2.41 


K 


5.6 


30.14 


0.35 


a 


OGLE-TR-132 


1.26 


1.34 


0.74 


F5 


3.5 


30.18 


0.1 


a 


OGLE-TR-56 


1.17 


1.32 


0.72 


GO 


4.2 


29.9 


0.2 


a 


TrES-1 


0.87 


0.82 


2.23 


KOV 


5.6 


30.04 


0.35 


a 


TrES-2 


0.98 


1.00 


1.38 


GOV 


4.2 


30.7 


0.2 


a 


TrES-3 


0.924 


0.81 


2.43 


G6V 


5 


30.38 


0.2 


a,b 


TrES-4 


1.384 


1.81 


0.33 


F8V 


3.8 


30.06 


0.1 


a 


WASP-1 


1.24 


1.38 


0.66 


F7V 


3.65 


30.12 


0.1 


a,c 


T 3 7 A CP* i~l 

WASP-2 


0.84 


0.83 


2.04 


KIV 


5.7 


30.00 


0.35 


a,c 


WASP-3 


1.24 


1.31 


0.78 


F7V 


3.65 


30.12 


0.1 


a,c 


WASP-4 


0.9 


1.15 


0.84 


G7V 


5 


30.38 


0.2 


a,c 


WASP-5 


0.972 


1.03 


1.27 


G4V 


4.8 


30.46 


0.2 


a,c 


WASP-6 


0.85 


0.9 


1.65 


G8 


5.1 


30.24 


0.2 


a,c 


WASP-7 


1.28 


1.24 


0.96 


F5V 


3.5 


30.18 


0.1 


a,c 


WASP-8 


0.87 


0.91 


1.63 


G6 


5 


30.28 


0.2 


a,c 


WASP-9 


1.1 


1.05 


1.34 


GO 


4.2 


30.6 


0.2 


a,c 


WASP-10 


0.71 


0.78 


2.09 


K5 


6.7 


29.7 


0.35 


a,c 


WASP-11 


0.74 


0.74 


2.58 


K3 


6.1 


29.94 


0.35 


a,c 


WASP-12 


1.1 


1.05 


1.34 


GO 


4.2 


30.6 


0.2 


a,d 


WASP-13 


1.3 


1.3 


0.84 


F9 


4 


29.98 


0.1 


a,c 


WASP-14 


1.319 


1.3 


0.85 


F5V 


3.5 


30.18 


0.1 


a,c 


WASP-15 


1.4 


1.4 


0.72 


F5 


3.5 


30.18 


0.1 


a,c 


XO-1 


1.00 


0.93 


1.77 


GIV 


4.3 


30.56 


0.2 


a 


XO-2 


0.98 


0.96 


1.54 


KOV 


5.6 


30.14 


0.35 


a 


XO-3 


1.213 


1.38 


0.66 


F5V 


3.5 


30.18 


0.1 


a 


XO-4 


1.32 


1.55 


0.5 


F5V 


3.5 


30.18 


0.1 


a 


XO-5 


1.00 


1.11 


1.03 


G8V 


5.1 


30.24 


0.2 


a 



Refer ences: htt£^^^www^exo 2la.net.eu| - retrieved 15/JUN/08. '' [O'Donovan et al.| ^mi\ 



et al. 



1 2008 I 



Hebb et al. 



1 2008 I 



Cameron 



value to the extremes does not affect our conclusions (see 
Section 3.3 1. 



For each of the exoplanets in our sample, we have cal- 
culated the total incident energy over the planet surface 
for the complete saturation interval of its host star. The 
stellar parameters required are listed in Table [l] and plane- 
tary properties in Table [2] The saturation duration is based 
only on the spectral type of the star. The X-ray saturation 
level used can be found by using the saturation level for 
stars of a given mass (in Table 3 of |Pizzolato et al.||2003"l ). 
Four stars in this work are slightly more massive than those 



studied by Pizzolato et al. (20031, and a saturated value 
of LJ"*/I/i,oi = 10~'^'^ is used, which is the average of the 
highest mass group in Pizzolato's work. 

A difference between our analysis and the population 



range of X-ray behaviour for stars of each type, whereas we 
assume average properties, since we do not know the indi- 
vidual histories of the objects in our study. Our estimates of 
X-ray irradiation would be biassed if the X-ray emission of 
stars were correlated with the presence of planets. Indeed, 



Kashyap et al. ( 2008 1 find evidence for excess X-ray emis- 



sion from stars with known planets. We choose not to modify 
our estimates by this factor because we are interested in the 
duration and strength of the saturated phase of X-ray emis- 
sion, whereas the measurements of Kasyap et al. apply to 
the weaker X-ray emission of much older stars. Nevertheless, 



if the results of Kashyap et al. ( 2008 1 were found to apply 



to stars of all ages, then we may have underestimated the 
X-ray irradiation of our sample by around half an order of 
magnitiude. 



study of Penz et al. ( 2008 1 is that Penz et al. included a 
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Table 2. Planetary parameters used in this work, includes calculated figures for (a), a-rocha PE and (^)- 
Figures for a, Rp and Mp are taken from The Exopla net Encyclop edia" unless otherwise specified. Densities 
may bo compared with that of Jupiter, 1.33 g cm~^ l | Allen|2000[ ) . 



Planet 


<a> 


^roche 


Mp 


Rp 


Pp 


\og{-PE) 




}roche 


Ref. 




(au) 


(au) 


(M,) 


(R.) 


(g cm"'"*) 


(ergs) 


(ergs r^/j) 


(ergs T-J^ ) 




HAT-P-lb 


0.055 


0.018 


0.53 


1.203 


0.38 


42.77 


42.33 


43.29 


a 


HAT-P-2b 


0.077 


0.006 


8.64 


0.952 


12.42 


45.29 


41.1 


43.21 


a 


HAT-P-3b 


0.039 


0.012 


0.599 


0.89 


1.05 


43 


42.44 


43.45 


a 


HAT-P-4b 


0.045 


0.018 


0.68 


1.27 


0.41 


42.96 


42.75 


43.51 


a 


HAT-P-5b 


0.041 


0.015 


1.06 


1.26 


0.66 


43.35 


42.62 


43.47 


a 


HAT-P-6b 


0.052 


0.017 


1.057 


1.33 


0.56 


43.32 


42.67 


43.65 


a 


HD 17156b 


0.184 


0.01 


3.08 


1.15 


2.51 


44.31 


40.65 


43.01 


a 


HD 189733b 


0.031 


0.012 


1.15 


1.154 


0.93 


43.46 


42.26 


43.09 


a 


HD 209458b 


0.045 


0.018 


0.69 


1.32 


0.37 


42.96 


42.79 


43.6 


a 


OGLE-TR-113b 


0.023 


0.011 


1.32 


1.09 


1.26 


43.6 


42.48 


43.13 


a 


OGLE-TR-132b 


0.031 


0.014 


1.14 


1.18 


0.86 


43.44 


42.21 


42.86 


a 


OGLE-TR-56b 


0.023 


0.015 


1.29 


1.3 


0.73 


43.51 


42.58 


42.94 


a 


TrES-lb 


0.039 


0.014 


0.61 


1.081 


0.6 


42.93 


42.6 


43.47 


a 


TrES-2b 


0.037 


0.011 


1.98 


1.22 


1.35 


43.91 


42.1 


43.12 


a 


TrES-3b 


0.023 


0.012 


1.92 


1.295 


1.1 


43.85 


42.25 


42.8 


a,b 


TrES-4b 


0.049 


0.023 


0.84 


1.674 


0.22 


43.02 


42.69 


43.33 


a,c 


WASP-lb 


0.038 


0.019 


0.89 


1.443 


0.37 


43.14 


42.93 


43.53 


a,c 


WASP-2b 


0.031 


0.012 


0.88 


1.038 


0.98 


43.27 


42.95 


43.75 


a,c 


WASP-3b 


0.032 


0.014 


1.76 


1.31 


0.97 


43.77 


42.21 


42.93 


a,c 


WASP-4b 


0.023 


0.016 


1.1215 


1.416 


0.49 


43.35 


42.32 


42.65 


a,c 


WASP-5b 


0.027 


0.011 


1.58 


1.09 


1.51 


43.76 


42.04 


42.81 


a,c 


WASP-6b 


0.045 


0.006 


1.3 


0.5 


12.8 


43.93 


40.8 


42.61 


a,c 


WASP-7b 


0.062 


0.012 


0.96 


0.915 


1.55 


43.4 


42.18 


43.61 


a,c 


WASP-8b 


0.082 


0.011 


2.35 


1.21 


1.64 


44.06 


41.07 


42.81 


a,c 


WASP-9b 


0.033 


0.011 


2.69 


1.22 


1.83 


44.17 


42.19 


43.14 


a,c 


WASP-lOb 


0.037 


0.009 


3.06 


1.29 


1.76 


44.26 


41.77 


42.96 


a,c 


WASP-llb 


0.043 


0.013 


0.58 


0.98 


0.76 


42.93 


42.24 


43.3 


a,c 


WASP- 12b 


0.0229 


0.021 


1.41 


1.79 


0.3 


43.44 


42.84 


42.92 


a,d 


WASP- 13b 


0.055 


0.017 


0.49 


1.06 


0.51 


42.75 


42.01 


43.04 


a,c 


WASP- 14b 


0.036 


0.008 


7.725 


1.259 


4.78 


45.07 


42.13 


43.4 


a,c 


WASP-15b 


0.048 


0.02 


0.56 


1.33 


0.42 


42.77 


42.52 


43.3 


a,c 


XO-lb 


0.049 


0.015 


0.9 


1.184 


0.67 


43.23 


42.59 


43.64 


a 


XO-2b 


0.037 


0.014 


0.57 


0.973 


0.77 


42.92 


42.57 


43.42 


a 


XO-3b 


0.047 


0.007 


11.79 


1.217 


8.08 


45.46 


41.87 


43.52 


a 


XO-4b 


0.056 


0.015 


1.72 


1.34 


0.88 


43.74 


41.8 


42.97 


a 


XO-5b 


0.051 


0.013 


1.15 


1.15 


0.93 


43.46 


41.41 


42.59 


a 



Refere n ces: °|http :/ /www .exoplanet. eu| - retrieved 15/JUN/08. 
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|0'Donovan et al.| | [2007l l. [Cameron et al.| 



3 RESULTS 



3.1 Past evaporation of known systems 

Figure [2] compares the potential energy of planets in our 
sample with the energy incident during the saturated X-ray 
emission phase of their parent star. As can be clearly seen, all 
of the systems considered in this work are able to survive the 
energy absorbed during the saturation period of their parent 
star. Of course, each of these systems has already weathered 
one saturation period, so Fig. |2] really tests whether these 
systems could survive a second saturation period. Since they 
all could survive another saturation period, it seems likely 
that they were only minimally affected by the first period. 
However, systems close to the dashed line may have lost 
around half of their mass during the saturation period they 
have aleady weathered. 



3.2 Evaporation at the Roche Hmit 

We know that all the planets in our sample have survived 
the saturation phase of their star. However most of these 
systems could have existed closer in, since most are well 
outside their Roche limit. In fact, there is a clear deficit 
of systems close to their Roche limits, as shown in Fig. [3j 
despite selection effects biasing surveys towards detection of 
close-in systems. The Roche limit sets the minimum orbital 
separation for a planet, within which tidal forces become too 
strong and pull the planet apart. 

The Roche limit for a fluid body (gas giant planets are 
assumed to be best represented by a fluid model rather than 
a rigid body), as given by |Aggarwal fc Oberbeck ( 19741, is 



aroche = 2.44a log 



(2) 
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Figure 2. A diagram of potential energy against the X-ray/EUV 
energy incident on the planet in one saturation period for the 
known transiting systems in their current configurations. The 
dashed line is an isocron representing one saturation period. 
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Figure 4. A diagram of potential energy against the X-ray/EUV 
energy incident on the planet in one saturation period for the 
known transiting systems, if they existed at their Roche limit. 
The dashed line is an isocron representing one saturation period. 
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Figure 3. A histogram showing the orbital separation of the 
sample of exoplanets from Table |2] scaled by their Roche limit. 
The majority of known systems can be moved approximately 2-3 
times closer to their parent star. 



where Uroche is the Roche limit, a is the semi major axis, and 
pp and p, are the density of the planet and star respectively. 

Re-plotting Figure [2] but moving the planets in to their 
Roche limit results in Figure |4] In this case ~50% of known 
systems would have absorbed enough energy to move their 
entire mass out to infinity. Therefore closer-orbiting analogs 
of the known planets would have suffered catastrophic evap- 
oration. The loss of such systems may account for the cut-off 
seen in Figure [T] since they would be situated to the left of 
their current positions, thereby filling the sparsely populated 
areas of these diagrams. Planets that would have survived 
at their Roche limit tend to have either high densities and 
surface gravities, or already exist close to their Roche limits. 



3.3 Effect on planet property distributions 

|Mazeh et all ( |2005| ) and |Southworth et al.| ( |2007[ ) have sug- 
gested that the apparent correlations found in mass-period 
and surface gravity-period plots for the known planetary sys- 
tems may be due to evaporation processes. We have shown 
that if planets like the known sample existed closer to their 



stars they would not have survived the saturated phase of 
stellar X-ray/EUV emission, and hence would appear as a 
void in these plots. 

If X-ray/EUV irradiation is important, then the un- 
derlying cut-off must be linear in the Mp/Rp versus 
plane, as the ratio of potential energy to irradiating energy 
per unit time defines the planet's characteristic lifetime. A 
plot of this form is shown in Figure [5] A sparsely populated 
region is seen in the lower left, similar to those in Figure 
[l] For a given class of star it is possible to plot a power 
law on this diagram, representing a destruction limit, be- 
low which planets would either have evaporated completely 
or left planetary cores that are not currently detectable in 
transit surveys. The functional form of this destruction limit 
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where Mp is the mass of the planet. Hp is the planetary 
radius, a is the semi major axis, LJ°* is the saturated 
X-ray/EUV luminosity of the star, Tsat is the saturation 
timescale of the spectral type and G is the gravitational 
constant. 

As can be seen in Figure [Sj the plotted lines correspond- 
ing to FGK stars track the edge of the missing region very 
well without any fitting. All systems exist above the destruc- 
tion line for their spectral type. 

Uncertainties and assumptions in our simple evapora- 



tion model (as discussed in Sections [2^ and 2.3 1 could move 
the destruction limits plotted in Fig.|5|, but they are unlikely 
to change the functional form, which matches the observed 
distribution well. 

The destruction limits can be raised by factors of 8, 
3 and 8.5 for F, G and K stars respectively before known 
planets (in their current state) fall into the danger zone. 
This suggests that our assumptions have tended to under- 
estimate the effects of evaporation. This is consistent with 
our expectations, since neglecting Roche geometry, planet 
radius evolution (Sect. [2^ , X-ray irradiation since the sat- 
urated phase, and the observed excess X-ray emission of 
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Figure 5. A diagram plotting mass squared over radius cubed 
versus the inverse of the semi major axis of the planetary or- 
bit squared. Planets that exist around G type stars are repre- 
sented with star symbols, K type stars are represented by straight 
crosses, and F type stars by diagonal crosses. The lines represent 
destruction limits, below which a system would have completely 
evaporated in the saturation time of its star. The dotted line 
represents this cutoff for GO stars, the line with large dashes rep- 
resents KO stars, and the line with small dashes represents F5 
stars. 



planet hosting stars (Sect. 2.31 are all likely to result in un- 



derestimates of evaporation rates. 

We note that uncertainties in saturation timescales 
alone cannot move the destruction limits up to the known 
planets. This would require timescales of 0.8, 0.6 and 3 
Gyr for F, G and K stars respectively, which are unreason- 
ably high values compared with existing X-ray observations 



(Sect. 2.31 



4 CONCLUSIONS 

We have shown that a simple model of exoplanet evapora- 
tion during the saturated phase of stellar activity is con- 
sistent with the observed population of transiting planets. 
The model also suggests that these same planets would not, 
in general, have survived at their Roche limit. We conclude 
that there is likely to have been a population of close-in 
planets that have been lost to evaporation. If this is the 
case there is no need for migration models to account for 
the observed lack of planets close to their Roche limits. 

This lost population can account for the correlations be- 
tween mass-period and surface gravity-period suggested by 
Mazeh et al.l (2005 ) 



and Southworth et al. 



( 2007[ ). Further 
to this we suggest that a cut-off in the M^/Rp versus 
plane is the underlying relation, and we show that evapora- 
tion in the saturated phase defines destruction limits that 
are a good match to the position of the lower edge of the 
observed distribution. 
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